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ABSTRACT 

Similarity solutions are found for the adiabatic collapse of density pertur- 
bations SM/M oc r~ s (s > 0) in a flat universe containing collisional gas only. 
The solutions are obtained for planar, cylindrical, and spherical perturbations 
with zero initial pressure. For adiabatic index 7 > 4/3, a shock develops at a 
fixed fraction of the current turnaround distance. Near the center of a spher- 
ical perturbations with 7 > 4/3 and s > 1/2, the gas is in quasi-hydrostatic 
equilibrium (pressure supported) and has an asymptotic power law density 
profile, p ~ r -3*/(»+i) t independent of 7. For s < 1/2, the profile depends on 
7, the pressure is finite, the temperature decreases inward, and gravity dom- 
inates pressure causing a continuous inward flow. Although for 1/2 < s < 2 
the temperature decreases at the center, the gas is pressure supported. The 
pressure is finite in cylindrical perturbations for s < 2(7 — l)/(&y — 4), and in 
planar perturbations for any s > 0. We also derive the asymptotic behaviour 
of the gas variables near the center in a universe dominated by collisionless 
matter. In such a universe, the gas in a spherical perturbation with s < 2 can- 
not be pressure supported and the temperature approaches a constant near 
the center. The solutions and the asymptotic behaviour are relevant for mod- 
elling the gas distribution in galaxy clusters and pancake-like superclusters, 
and determining the structure of haloes of self-interacting dark matter with 
large interaction cross section. 

Key words: cosmology: theory - gravitation - dark matter -baryons- in- 
tergalactic medium 



1 INTRODUCTION 
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On scales larger than a few megaparsecs, pressure forces 
in the baryonic matter in the universe are negligible, so 
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the evolution of dark and baryonic matter is mainly 
determined by gravity. On small scales pressure be- 
comes important and may segregate between the evo- 
lution of baryonic and dark matter. Pressure forces, 
cooling of gas, and star formation feedback, are key in- 
gredients in galaxy formation. These ingredients com- 
bine to cause differences between the distributions of 
galaxies and dark matter (biasing), even on large scales 
where these effects are not directly important (e.g., 
Kaiser 1984, Dekel & Rees 1987, Kauffmann, Nusser 
& Steinmetz 1997, Benson et. al. 2000). On scales 
smaller than the Jeans length of the photo-heated in- 
tergalactic medium (IGM), pressure forces dominate 
gravity and can prevent the collapse of gas into dark 
haloes below a certain mass threshold. For haloes mas- 
sive enough the temperature of the IGM can be ne- 
glected and the gas falls into the halo. The mean free 
path for collisions between gas particles inside a halo is 
» (200/5 c ) (l0- 15 /a) (1-Mr 3 (0.1/n b h 2 ) 1.6pc, where 
5 C is the overdensity inside the virial radius and a is a 
typical cross section for collisions in units of cm 2 . This 
is smaller than the virial radius of a typical halo by a 
few orders of magnitude. Therefore, on its infall into 
the halo, the gas is likely to form shocks and trans- 
form its kinetic energy into heat. The hot dense gas can 
then cool to form stars which explode and inject energy 
into the halo gas. Detailed study of these processes un- 
der general conditions is not feasible. One can aim at 
a global parameterization based on general physical re- 
quirements which match observational data (Kauffman, 
White, Guiderdoni 1993, Somerville & Primack 1999, 
Cole et. al. 1994). Another route would be to study 
special aspects which can be treated by either numeri- 
cal or analytical methods. Here we focus on the collapse 
of the baryonic gas in an Einstein-De Sitter universe, ig- 
noring the gas initial temperature, cooling and heating 
processes. We assume that the collapse initiates from 
a symmetric scale free density peak, and that the ve- 
locity of each shell in the peak is taken to match the 
general expansion of the universe. The energy of each 
shell is negative and it will expand up to a maximum 



distance before it starts falling towards the center of 
the perturbation. The maximum distance is termed the 
turnaround radius. Shell crossing is not allowed and the 
collapse can proceed in two distinct ways, either a shock 
wave forms, or shells accumulate at the center. Which 
of these possibilities actually occurs, depends on the 
physical conditions at the center. If the velocity van- 
ishes at the center than a shock wave forms. If on the 
other hand physical conditions allow a non vanishing 
velocity at the center then the shells accumulate at the 
center (Bertschinger 1985). Inner boundary conditions 
can be arranged so that a shock is accompanied by 
the accumulation of central mass, but a proper stabil- 
ity analysis is needed to determine whether or not this 
is possible (Bertschinger 1985). In this paper we focus 
on shocked collapses without the formation of a central 
mass. Since the initial gas pressure is negligible, the 
collapse eventually develops in a self-similar way where 
the only relevant scale at any time is the radius of the 
shell at maximum expansion. Bertschinger (1985), and 
Forcada-Miro & White (1997) have studied similarity 
solutions in spherically symmetric perturbations with 
initial relative mass excess SM/M oc r -3 , and r~ 2 , re- 
spectively. Here we derive similarity solutions in planar, 
cylindrical, and spherical geometries, for the collapse of 
a perturbation with 5M/M oc r~ s for any s > and 
adiabatic index 7 > 4/3. 

In section |^ we write the equations of motion 
for symmetric perturbations in planar, cylindrical, and 
spherical geometry. In section ^ we discuss the asymp- 
totic behaviour of the fluid variables near the center, 
in the case of shocked collapse. In section ^ we present 
results of numerical integrations of the equations. In 
section ^| we derive asymptotic behaviour of the fluid 
variables in a universe dominated by collisionless dark 
matter. In section ^ we conclude with a discussion of the 
results and their potential astrophysical consequences. 
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2 THE EQUATIONS 

We write the Newtonian equations of motion governing 
the adiabatic collapse of symmetric perturbations in a 
collisional fluid (gas) of adiabatic index 7 and zero ini- 
tial pressure. Except section ^, we restrict the analysis 
here and throughout to the collapse in a flat universe 
containing collisional gas only. The initial gas pressure 
is zero, so the expansion scale factor of the universe 
is ait) oc t 2/3 , the Hubble function is H(t) = 2/(3t), 
and the background density is p c = 3H 2 /(SttG) = 
l/(67rGi 2 ). 

Denote by r and v = dr/dt the physical position 
and velocity of a gas shell, where r — is the sym- 
metry center of the perturbation. Further, let p(r, t) 
and p(r,t) be the gas density and pressure at r. As 
in Fillmore & Goldreich (1984) define the mass within 
a distance r from the symmetry center by m(r,t) = 
J x"~ 1 p(x, t)dx, where n = 1,2, and 3 refer, respec- 
tively, to planar, cylindrical, and spherical perturba- 
tions. The mass within a fixed shell varies with time 
like m ~ t- 2<3 - n)/3 , because of the Hubble expansion 
along 3 — n of the axes. In this notation, the equations 
of motion are, the continuity equation, 



to zero, ti, as 
8M 



2(3-n) . 



d{pt~^— ) _ .aa^si i_ nn ,_„_i. 



At 
Euler, 

dv 2 3 — n r 



■ts—pr L - n d r (r n - L v) 



dt 9 n t 2 
adiabatic condition, 



s (pp-) = o 



and the relation, 



o n— 1 

a r m = r p . 



d r p 4irGm 

p r n-l 



(1) 



(2) 



(3) 



(4) 



These equations are satisfied everywhere except at 
the shock where the fluid variables are described by 
jump conditions obtained from mass, momentum, and 
energy conservation. The initial conditions leading to 
self-similar collapse are specified at an early time close 



SM _ /M~ s 
M \r ) ' 

2 

v(r,ti) = —r, 
p{r,U) = 0, 



(5) 

(6) 
(7) 



where 8M/M is the mean density contrast interior to r, 
and s > 0. For cosmological initial conditions the initial 
density contrast must be tiny, so we will be interested in 
the solution in the region r 2> Tq. A perturbation with 
s > n can be realized by placing a high narrow posi- 
tive density peak at the center (<C ro) of a symmetric 
void with local density contrast ~ (— r~ 3 ). The condi- 
tion (|^) means that a gas shell at r moves initially with 
the general universal expansion. This condition can be 
relaxed to allow for a non vanishing initial zero pecu- 
liar velocity according to late time linear theory (e.g., 
Peebles 1980). However, this does not affect the details 
of the collapse (Peebles 1980, Bertschinger 1985), so we 
use (^) which is commonly adopted in the literature. 
Bertschinger (1985) and White & Forcada (1997), re- 
spectively, considered the collapse of spherical pertur- 
bations with s — 3, and s — 2. 

The equations of motion (Q-^) together with the 
initial conditions (^[-Q) are insufficient to completely de- 
termine the evolution of the perturbation. Still missing 
is an inner boundary condition specifying the velocity 
and mass at r = 0, for t > ti. For a shock to develop 
without the accumulation of a central mass (a black 
hole for n — 3) we must have v(r = 0, t > ti) = and 
m(r = 0, t > U) — 0. Relaxing the condition V(0) = 
leads to a non-vanishing central mass with or without 
the presence of a shock. 

In a critical density universe (fi = 1) the only 
length scale relevant to the collapse is the scale of non- 
linearity. At any time, t, this scale can be defined as the 
distance of the shell at the maximum expansion, i.e., 
the shell with v — (Gunn 1977, Fillmore & Goldreich 
1984, Bertschinger 1985) . This radius is termed the cur- 
rent turnaround radius, r ta (t). Starting from tiny initial 
density contrast, the mean overdensity (density in units 
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of p c ) interior to rt a (t) is a fixed number independent 
of time. For time t S> ti, when shells with r > ro reach 
their turnaround, the collapse develops a self-similar be- 
haviour that depends on r and t through the combina- 
tion A = rjtta- The turnaround radius rt a {i) is given 
by (e.g., Fillmore & Goldreich 1984), 



r ta = r - 



c x ft_\ a 



c*; 



2 s + 1 
a = 

3 s 



(8) 



where, 



3tt\ 2 / 3 



(+=^,0.741,1 ; C t = jj, 1.386, (!M 



for n = 1, 2, and 3, respectively. The turnaround radius 
grows faster than the scale factor a ~ t 2//3 . This is be- 
cause the mass, ~ p c »~ta(i) 3 , interior to r ta must grow 
with time while the mass, ~ p c a 3 (t), inside a fixed shell 
in a homogeneous universe is constant. For s < 2 the 
turnaround radius grows faster than t reaching the hori- 
zon scale in finite time. When this happens relativistic 
description must be used and rt a ceases to be the only 
scale in the problem (Fillmore & Goldreich 1984). 

The equations can be cast into a non-dimensional 
form using the scaled variables V(A), D(\), -P(A), and 
M(A) defined by (Bertschinger 1985), 
rta , 



v(r,t) 



t 



-V(X) 



p(r,t) = p c D(X) 



p(r,t) 



Pc 



(^) 2 P(A) 



m(r,t) = -p c r" a M(X) 



(10) 

(11) 
(12) 

(13) 



Expressed in terms of these variables, the equations (Q- 
^ become, respectively, 



(V-a\)D'+ (2—±v + V' -?f)D = 



(a - 1) V+CV - aX) V A = — — — - — 

v ' x ' 9 n D 9A"- 



M' = 3A n_1 L> 



(14) 



(16) 
(17) 



We will mainly be concerned with solutions for 
shocked collapse with vanishing mass at the center. The 
inner boundary condition appropriate for this collapse 
are vanishing mass and velocity at A = 0, i.e., 



V(0) = and M(0) = , 



(18) 



Self-similarity implies that the shock appears at fixed 
X = X s — r s /r ta , so the physical radius of the shock r s oc 
t a and its non-dimensional speed is (r ta /t)~ 1 (dr B /dt) = 
aX s . At the surface of the shock the fluid variables 
satisfy the jump conditions obtained from mass, mo- 
mentum, and energy conservation. In terms of the non- 
dimensional fluid variables, the jump conditions appro- 
priate for an adiabatic shock are, 
7-1 



V + = aX s + 



7+1 



(V- - aX s ) , 



7 + 1 
M + = M~ , 



7-1 
1 D~(y--aX s f 



(19) 
(20) 

(21) 
(22) 



where the superscripts of the minus and plus signs re- 
fer to pre- and post-shock quantities. In employing en- 
ergy conservation we have taken QT^rn as l ne non- 
dimensional internal energy per unit mass. 

In section ^|we will find numerical solutions satisfy- 
ing the requirements for shocked collapse without a cen- 
tral mass. Except spherical perturbations with 7 = 4/3 
only one value A s can yield solutions satisfying these re- 
quirements. Spherical perturbations with 7 = 4/3 allow 
a range of values for A s . Before presenting the numerical 
solutions we derive in the next section the asymptotic 
behaviour of the fluid variables near the center, and two 
integrals of motion which will be used as a check on the 
numerical solutions. 



where the prime symbol denotes derivatives with re- 
spect to A. 
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3 INTEGRALS OF MOTION AND 

ASYMPTOTIC BEHAVIOUR NEAR THE 
CENTER 



Solutions to (J14| — |l 7|) with the appropriate jump and 
boundary conditions for all A will be found by numerical 
integration. We present here an analytic treatment of 
the equations to derive the asymptotic behaviour of the 
fluid variables near r = 0, and two integrals of motion 
(e.g., Bertschinger 1983, 1985). We restrict the analysis 
shocked collapses satisfying the inner boundary condi- 
tion @. All fluid variables can be expressed in terms 
of an auxiliary function K(X) = exp 

Jo vw-tt* as fol " 

lows. 



K 



V(X) = ^7+«A 

D 
P(A) 



■ o 



M(X) 



(23) 

X l-n K , K -l + n (2/3- a ) ^ (24) 
^(l-«)7^-4+(2n/3-3)7-a(2+n7) j^t~j /ggs 

(26) 



gJ(n(2/3— a) 



D n(3a - 2) 

where Do and Po are constants. The fluid variables 
in (p3j-p6|) satisfy the non-dimensional equations (|l4|), 
(^), and ( |l7j ) for any functional form of K. The func- 
tion K is then specified by only one equation, the non- 
dimensional Euler equation (jt5|). 

Two integrals of motion can immediately be found 
from (]23|-p6]). These are the mass and entropy integrals 
of motion (Bertschinger 1983, 1985), 
9 



M 

PD-"'M C 



n(2 - 3a) 
const ; 



D(V -aX)X" 



(27) 



(28) 



n 2 — 3a 

where all fluid variable are evaluated at any A inside 
the shock. Since £ < 0, the entropy integral of motion 
means that the entropy tx ln(PD~ 7 ) is an increasing 
function of the mass and hence A. 

The auxiliary function greatly simplifies the deriva- 
tion of the asymptotic behaviour of the non-dimensional 
fluid variables near A = 0. Since V(0) = and M(0) = 
0, equations (^3|) and (|2(]) imply that, to first order, 
K(X) must approach 
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K(X) = X 1/(v °- a) (29) 

as A — > 0, where Vo is an arbitrary constant. Substitut- 
ing this expression for K in ( p3| - ^) , yields 



V = V X 
M = 



9A> A ^ 



D 
P 



n{3a - 2)' 
Do A , 



a - V 
Po 



-A' ; 



(30) 
(31) 

(32) 

(33) 



(a - Vo)'' 

and the asymptotic exponents are expressed in terms of 
the coefficient Vo as 
. _ n(2 - 3V ) 



n 



3(V -a) ' 

4-2a-2 7 +(§- V )n-y 



(34) 
(35) 



Vo — a 

These relations have been obtained without using the 
Euler equation (|l5|). In order to determine the expo- 
nents uniquely we use the Euler equation which adds 
the following constraints, 



77 = 2(5 + 1), for r? < , 
2(8 + 1) > 0, for t? = 



(36) 



A value 77 > gives Po < and so, by (Eq), a negative 
entropy integral of motion. But the jump condition ( pl[ ) 
gives positive pressure, P + , just behind the shock and 
since the density and mass are also positive, the entropy 
integral (|i^) must be positive. So we rule out 77 > 0. 

If the solution to @, @, and @ is rj < then 
the Euler equation also provides the following constraint 
on the coefficients Po and Do, 



Po = 



2P 2 a(l/ -a) 
n?7(3a — 2) 



7-1 



(37) 



Table 1 lists the values of Vo, 77, and 5 in all cases. The 
dashed curves in figure |l] are a graphical representa- 
tion of the density exponent 5 versus s for 7 = 5/3. 
In planar geometry, n = 1, the only possible solution to 
equations is 77 = and the pressure is finite everywhere. 
In cylindrical geometry, n = 2, we have 77 < only for 
a < 3(7-1) ' so ^ mus t be zero for other values of a. In 
spherical geometry, n — 3, if a < 2 and 7 > 4/3 then 
77 = 2(1 — — ), and Vo = meaning that V(X) oc X v , 
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where v > 1. A second order expansion gives v = 1 -+- 



(4 - 5a) ± 7 - 1/2 ^80(7 - 8) + I67 + a 2 (32 + 7) 
where only one of the roots is v > 1. In the limit of ei- 
ther a — > 2 or 7 — * 4/3, we have v — > 1. 



2a 



For n = 3 and 7 = 4/3 the relations (|3J), (|35j), and 
(^) allow multiple solutions for Vo and consequently 
for A s . A second order calculation gives the upper limit 
Vo < 4 — 5a. Solutions exist for any positive A s smaller 
than a maximal value which corresponds to the upper 
limit on Vo. This means that for 7 = 4/3 a shocked 
collapse cannot be accompanied by the presence of a 
non- vanishing mass at the center. In the table we list the 
asymptotic constants corresponding to Vo = 4 — 5a, i.e., 
the maximal value of A s for which a shocked collapse 
occurs. Bertschinger (1984) does not mention that there 
are solutions for shocked without a central mass for a 
range of A 3 . His numerical solution with 7 = 4/3 seems 
to correspond to the maximal A s and yields 77 = —3.2 
and 5 = —2.6, instead of —3 and —2.5 as listed in table 
1. 

We now examine how the dimensional density, 
p(r,t) — p c D, varies with time near the center. Using 
p c ~ t~ 2 , the first order expression D ~ (r/rt a ) S , and 
r ta ~ t a , we find p(r,t) ~ r 6 /t 2+aS . When 2 + 5a = 
the density is constant with time. In spherical geome- 
try a time independent density is equivalent to Vo = 0. 
Because of the expansion in the n — 1 directions, a van- 
ishing Vo in planar and cylindrical geometries leads to 
p ~ f -2-a<5 ^ t 2(n-i)/3_ In all geome tries a vanishing Vo 

indicates that, to first order in the asymptotic expan- 
sion, the gas is pressure supported in hydrostatic equi- 
librium. According to the table, spherical perturbations 
have Vo = only for s > 1/2. Note that for 2 > s > 1/2 
the asymptotic density is constant even though the tem- 
perature ~ P/D ~ r v ~ s decreases inward. For s < 1/2, 
the density increases with time. According to table 1, 
planar and cylindrical perturbations have Vo = only 
at s = 1/(2 — 7) and s = 2/(4 — 7), respectively. 



4 NUMERICAL INTEGRATION 

We present results of the numerical integration of the 
non-dimensional equations of motion. The numerical so- 
lutions shown here describe shocked collapses without a 
central mass. Outside the shock the fluid variables are 
given from the solution for collapse with zero pressure 
(Zel'dovich 1970, Peebles 1980, Fillmore & Goldreich 
1984). 

The shock position A s is unknown a priori. We have 
to find its value such that the fluid variables satisfy the 
equations of motion, (p^-|l^), the jump conditions, ([l9|- 
p2]), and the inner boundary condition, ([l8|). Assum- 
ing that the pre-shock variables are given from the zero 
pressure solution, the value of A 3 can be found as follows 
(Bertschinger 1985, Forcada-Miro & White 1997). For 
an assumed value for A s , we obtain the post-shock vari- 
ables using the jump conditions. We then integrate the 
equations of motion from A s inward to A = 0, and tune 
A s so that the solution gives V = and M = at A = 0. 
In all numerical solutions we find that if M(0) = for 
a given A 3 then V(0) = 0, and vice verse. 

The zero pressure solutions in planar and spheri- 
cal geometries are known analytically ( Zel'dovich 1970, 
Fillmore & Goldreich 1984) so the pre-shock fluid vari- 
able can be found directly, for an assumed A s . How- 
ever in cylindrical geometry an analytic solution is not 
available and we numerically integrate the equations 
with zero pressure to obtain the pre-shock quantities. 
In practice we use numerical integration also in pla- 
nar and spherical geometries. It is convenient to inte- 
grate the zero pressure equations from A = 1, i.e., the 
turnaround radius, to A s . At A = 1 the fluid variables 
are 

1/(1) = , P(1) = 0, (38) 

where the expressions for D(l) and M (1) were obtained 
from the zero pressure solution (e.g., Fillmore & Gol- 
dreich 1984). 
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The sequence of figures (Q-q) shows the numeri- 
cal solutions for the over density D, the pressure P, 
the velocity V, and the thermal energy U = o7~n m 
spherical, cylindrical and planar geometries for several 
values of s, all as a function of A. All curves are ob- 
tained from the solutions with 7 = 5/3. The solid curve 
in each plot corresponds to s = n. The sudden change in 
the fluid variables indicate the location of the shock. All 
numerical solutions satisfy the integrals of motion ^t\) 
and ( ps| ) up to the numerical accuracy and agree with 
asymptotic behaviour of the previous section. The loga- 
rithmic slopes of all curves match the corresponding val- 
ues listed in table 1. Spherical perturbations with s = 3 
and s = 2 were, respectively, analyzed by Bertschinger 
(1984) and Forcada-Miro & White (1997). The agree- 
ment between their solutions with 7 = 5/3 and ours is 
excellent. In all geometries the density and pressure in 
the shock are higher for larger s. Figures ^ of the ve- 
locity and figure [| of the internal energy U = p7~n 
demonstrate that particles are decelerated at the shock 
converting most of their kinetic energy into heat. The 
velocity near the center of cylindrical perturbations can 
be positive inside the shock, hence the linear vertical 
scale in velocity plot in this case. 

In figure ^ we plot the location of the shock X a as 
a function of s for several values of 7. The shock loca- 
tion always increases with s and agrees with the values 
obtained by Bertschinger (1985) and Forcada-Miro & 
White (1997) for spherical perturbations with s = 3 
and s = 2, respectively. A special case is spherical col- 
lapse with 7 = 4/3. Here a solution for shocked collapse 
without a central point mass is possible for any A s less 
than a maximal value A c . In this case we plot the maxi- 
mal value A c . For s ^> 1, A s can exceed unity increasing 
to a finite value as s — » 00. 

The variation of the fluid variables attached to a 
given fluid elements are also of interest. Using the ve- 
locity obtained from the self-similar solutions we can 
calculate the trajectory r(t) of a fluid element (parti- 
cle) as a function of time. The fluid variable associated 
with the particle at any time can then be obtained by 
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interpolating the solutions at the particle's position at 
that time. The three panels from top to bottom in figure 
show, respectively, the trajectory, density, and pres- 
sure of a particle as obtained from the solutions with 
7 = 5/3 for three values of s. The time axis in all pan- 
els is t/tta where tta is the time at which the particle 
reached its maximum expansion (the turnaround time) . 
The particle position, r, density p, have been scaled by 
their values at t ta , r' ta = r(tt a ), and p = p(tta), respec- 
tively. The pressure, p, has been scaled by its value at 
the shock p s . With this scaling the curves in the figures 
become valid for all the particles. The particle trajec- 
tory for s — 3 (solid line, top panel) agrees with the 
corresponding curve in figure 4 of Bertschinger (1985). 
As expected from the asymptotic solution the particles 
for s = 3 and s = 10 tend to settle at a physical dis- 
tance which is a fixed fraction of their turnaround radii. 
The decay of the trajectory for s = 1/4 at late time can 
be evaluated using the asymptotic expansion. Near the 
origin V(A) = VoA, so dr/dt = {r t a/t)Vo\ = —r/t and 
r ~ t v °. According to table 1 Vo = -8/15 for s = 1/4 
and 7 = 5/3 so r decays like t~ s ^ 15 . The density (mid- 
dle panel) and pressure (bottom) curves for s = 3 and 
s — 10 flatten at late times, consistent with the settling 
of the particles to a constant r. 

5 ASYMPTOTIC BEHAVIOUR IN A 
UNIVERSE DOMINATED BY 
COLLISIONLESS MATTER 

So far we have considered similarity solutions for col- 
lapse involving gas only. The collapse of scale free sym- 
metric perturbations in a an Einstein-De Sitter universe 
containing a mixture of gas and collisionless matter is 
also self-similar. Similarity solutions for mixed collapse 
are beyond the scope of the present paper. Here we only 
obtain the asymptotic behaviour of the gas variables in 
a universe dominated by collisionless matter. For a self- 
similar collapse to develop, the two matter components 
must start with the same initial conditions with zero 
initial gas pressure. Shells of gas and collisional matter 
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then move together until they reach either a shock in 
the gas or the region of shell crossing in the collisionless 
component. The evolved density profiles of both com- 
ponents depend on r and t through A. Although the 
global gas mass fraction is negligible, the gas can be 
gravitationally dominant at the center of the collapse 
if it has a steeper density profile than the collisionless 
matter. For the purpose of deriving the asymptotic ex- 
ponents we proceed assuming that the gravity of the 
gas is negligible everywhere and check the consistency 
of this assumption according to the results. So we re- 
place the gas mass M (A) in the non-dimensional Euler 
equation ( |l5| ) by the collisionless matter mass M(A). 
Writing M ~ \ S+n near A = 0, Fillmore & Goldreich 
(1984) find that 6 is 



n = 1 



8 = 



n = 2 : 



n — 3 : 



s + 1 
5 = -l 

5 = -2 , for s < 2 

for s > 2 



(40) 
(41) 

(42) 

s+1 ' v ; 

Table 2 summarizes the values of the asymptotic con- 
stants obtained by substituting the asymptotic be- 
haviour of the collisionless mass in the non-dimensional 
Euler equation. The dotted and dashed curves in figure 
[j] represent the density exponent vs s computed from 
table 1 and table 2, respectively, for 7 = 5/3. The solid 
curve is the collisionless matter exponent 8 versus s. 

In spherical perturbations we see from the tables 
and figure [l] that the asymptotic gas density profile in 
a universe with collisionless matter is steeper than that 
containing gas only. For s > 2 the density asymptotic 
exponents of the gas and collisionless matter are iden- 
tical. For s < 2 the collisionless matter density profile 
is steeper. This is consistent with neglecting the grav- 
ity of the gas near the center. As 7 — > 4/3 gas density 
exponent approaches —2 for s < 2 and — 3s/(l + s) oth- 
erwise, so the gas has the same asymptotic profile as 
the collisionless matter for all s. 



In cylindrical perturbations with s large enough, 
the gas has a steeper density profile than the collision- 
less matter. This means that virial motions in the col- 
lisionless matter are more effective at balancing gravity 
than the gas pressure force. So for large s the derivation 
of the asymptotic exponent in section ^| is more suitable. 
A steeper gas profile also occurs in planar perturbations 
with s < 1/(3 — 7). However the asymptotic constants 
obtained here and in section ^ are identical in planar 
geometry. 

In spherical perturbations with s > 2, the den- 
sity p{r,t) is constant with time and the temperature 
~ p/p diverges like r 1 + s as r — * 0. To first order in the 
asymptotic expansion where the gas can be in hydro- 
static equilibrium in the dominant potential well of the 
collisionless matter. For s < 2 the density increases with 
time and the temperature is constant near the center. 
For comparison, in the absence of collisionless matter, 
for 2 > s > 1/2 the temperature decreases inward but 
the gas is in hydrostatic equilibrium to first order in the 
asymptotic expansion. 



6 DISCUSSION 

The similarity solutions are found for collapse in a flat 
universe with matter density parameter Q = 1. Be- 
cause of Birchhoff's theorem, the solutions for spher- 
ical collapse are valid in an open universe if the current 
turnaround radius is well inside the spherical region in- 
terior to which the perturbation is bound. The state- 
ment is incorrect for planar and cylindrical perturba- 
tions because of the explicit appearance of cosmology 
dependent terms in the equations of motion (|l|-Q), like 
t 2(3-n)/3 in the cont i nmty equation (jl|). 

The solutions are appropriate for the adiabatic col- 
lapse of perturbations with deep gravitational poten- 
tial so that the initial thermal energy of the gas can 
be ignored. Such perturbations are probably the seeds 
for massive galaxies, galaxy groups and clusters. In the 
intergalactic medium (IGM) most of the gas is continu- 
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ously photo-heated and is of moderate density. There is 
considerable interest in analytic modelling of the IGM 
in current methods for extracting cosmological informa- 
tion from the Lyman forest (Croft et. al. 1998, Nusser & 
Haehnelt 1999, 2000). So far these methods have heav- 
ily relied on linear analysis (e.g., Bi, Borner, & Chu 
1992, Gnedin & Hui 1998, Nusser 2000) and hydrody- 
namical simulations (e.g., Petitjean et. al. 1995, Theuns 
et. al. 1999). Analytic treatment of the IGM beyond the 
linear regime is exceedingly complicated. Consider a sit- 
uation in which photo-heating establishes the relation 
p — kp 1 in the IGM, where k and 7 depend non-trivially 
on time (e.g., Theuns et. al. 1999). The pressure in this 
case introduces a length scale fc 1 ^ 2 G^ 1 fi/ 2 t 2 ' (e.g., 
Sedov 1959). If we take constant k and 7, this length 
scale varies with time like rt a only in the special case of 
7 = 4/3 and infinite s. So physically interesting situa- 
tions including initial pressure in which the collapse is 
self-similar do not exist. 

Spherical Perturbations with s > 2 when 7 > 4/3, 
and with any s > when 7 = 4/3 deserve special at- 
tention. In the corresponding solutions for shocked col- 
lapse without a black hole at the center, the quantity 
V c 2 = 2Gm/r diverges towards the center, and so there 
is a point r — r g at which V c 2 = c 2 , where c is the speed 
of light. According to general relativity this implies the 
presence of a black hole at the center, invalidating the 
assumption of no central mass, made in deriving the 
solutions. In particular, the condition v(r = 0) = 0, 
necessary for shocked collapse, is incompatible with the 
presence of a central black hole. Near r g , however, radi- 
ation pressure and angular momentum can prevent the 
formation of a black hole. Should this occur, we expect 
our solution for shocked accretion to be valid away from 
the central region. 

The evolved gas variables in the symmetric self- 
similar collapse contain full information on the initial 
perturbation. So the system retains memory of the ini- 
tial conditions, even in the highly nonlinear regime. On 
the other hand, a collapsing system of collisionless mat- 
ter can develop density profiles which do not depend on 
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the initial shape of the perturbation. For example, ac- 
cording to the solutions of Fillmore & Goldreich (1984), 
a spherical density perturbation develops into r~ 2 for 
s < 2, and a cylindrical perturbation into r for all 
s. Haloes identified in cosmological simulations of col- 
lisionless particles with generic initial conditions, also 
tend to have density profiles independent of the spec- 
trum of the initial fluctuations (Navarro, Frenk & White 
1997). 

Our results are relevant for describing the gas dis- 
tribution in various physical systems such as the cores 
of clusters or pancake-like superclusters. Over a limited 
range of scales, the index s can be related to the index, 
/, of the three dimensional power spectrum, p(k) ~ k l , 
of the linear density fluctuations. If the initial density 
field is gaussian with a scale free power spectrum then 
the properties of the nonlinear field depend only on one 
scale. This is the nonlinear scale, R n i, defined as the 
scale on which the rms value of density fluctuations is 

rn 2(1 + 5) 

unity. This scale grows with time liken R n i ~ t 3 ( l + 3 "> . By 

2(3+1) 

matching the time dependence of R n i and r ta ~ t 3s 
we identify s = (I + 3) /2. So the collapse of gas into 
galaxies and clusters can be, respectively, modeled by 
our solutions for s~0.4 and 0.65, where we have taken 
i~— 2.2 and —1.7 assuming a Cold Dark Matter power 
spectrum. Taking i~— 1 for collapse on a pancake-like 
large scale superclusters gives s ~ 1. Another way to 
relate s and I is to identify symmetric perturbations 
with local maxima in the linear density field (Hoffman 
& Shaham 1985). The shape of high density peaks in a 
gaussian field varies with r like the two-point correla- 
tion function, ~ r _ ( i+3 ) (e.g., Bardeen et. al. 1986). So, 
at least in the limit of high peaks, s = I + 3. On galaxy 
and cluster scales the relation s — 1 + 3, respectively, 
gives s = 0.8 and 1.3 in contrast to s = 0.4 and 0.65 
obtained from s = (I + 3)/2. Note however that gas in 

* R n l does not involve the dimension n because I refers to 
the three dimensional p(k) so the rms value on a scale R is 
R( l + 3 )/ 2 independent of n. 
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galaxies tends to settle into disks and so our solutions 
are less relevant than they are for cluster size objects. 

In spherical geometry the asymptotic behaviour 
shows that the gas cannot be pressure supported if 
s < 1/2, and s < 2 for collapse with, and without col- 
lisionless matter, respectively. Estimates of the masses 
of rich galaxy clusters from X-ray observations of the 
intracluster gas rely on hydrostatic equilibrium (e.g., 
Fabian 1994). If on cluster scales s ~ 0.7— ~ 1.3, then 
the asymptotic behaviour implies that the cluster gas 
may not be in hydrostatic equilibrium. How large is 
the error introduced in the mass estimates by assuming 
hydrostatic equilibrium? The following argument shows 
that this error is negligible. Hydrostatic equilibrium cal- 
culations neglect the term G~ 1 r 2 dv/At in the mass es- 
timate. Using the asymptotic expansion one finds that 
neglecting this term amounts to a relative mass error of 
~ {2/iv 2 )(t/tta) iVo ~ 2 where U a is the turnaround time 
of the shell present at r at the current time t. Shells in 
the inner regions have passed their maximum expansion 
a few dynamical times ago. Therefore t S> t ta and since 
Vo < we conclude that the error is negligible. 

The solutions are related to modelling the structure 
of haloes made of self interacting dark matter (SIDM) 
(Spergel & Steinhardt 1999) with large interaction cross 
section. On scales of massive galaxies and clusters, our 
results predict final density profile p~r -1,2 to ~r _1 ' 7 . 
These profiles are consistent with the results obtained 
by Moore et. al. in their simulations of SIDM with large 
cross section. 
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Table 1. Asymptotic constants, Vo, <5, and r\ for collapse 
of collisional gas only. Spherical perturbations with 7 = 4/3 
allow a range of A s . Listed are the values corresponding to 
the maximal A s . 

n = 3, 7 > 4/3 







V 




S 


s < 1/2 




4(2s-l) 





-3(8+s(37-4)) 




9s 7 


2+3 7 +s(37-4) 


s > 1/2 







r\ l — 2s 


3s 
s + 1 


n = 3, 7 = 4/3 






2(a-S) 
3s 


-3 


-2.5 


n = 2, 7 > 4/3 


s < 2 <?~ 


i) 


s(4 — f )-2 





-2(2+s(3 7 -4)) 


— 37- 


4 


3s7 


2+2 7 +s(3 7 -4) 


^ 2(7- 


1) 


2-7 


4(l-7) + 2s(37-4) 


2s(4-37) 


37- 


4 


3(7-1) 


2(l-7) + s(4-3 7 ) 


2(7-l) + s(37-4) 


n = 1, 7 > 4/3 






-4(l + s(7-2)) 





s(4-3 7 )-2 








2+7+s(3 7 -4) 


Table 2. 


The 


gas asymptotic constants, Vb,5, 77 for mixed 


collapse with dominant collisionless matter. 










n = 3, 7 > 4/3 








V 


»7 


5 


s < 2 




2(s-2) 


2+ s (3 7 -4) 


2+s(3 7 -4) 




9s(7-l) 


(l+4s)/3-7(l+s) 


(l+4s)/3-7(l+s) 


s > 2 









3s 
s + 1 


n = 2, 7 > 4/3 


« < 2( r 


i) 


B (4 — f )-2 





-2(2+s(37-4) 


— 37- 


4 


3s7 


2+2 7 +s(3 7 -4) 


. > 2 ir 


l) 


2(l+ 3 ( T -3)) 


2( 7 -l)-s(37-4) 


-2+2s(4-37) 


37- 


4 


3s(l-2 7 ) 


27+s(37-4) 


27+s(37-4) 


n = 1, 7 > 4/3 






-4(l + s(7-2)) 





s(4-37)-2 






3s7 


2+7+s(3 7 -4) 





-1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 



-0.5 : 

-1 K 
-1.5 
-2 



n=3 



1 1 



Figure 1. The density asymptotic exponent versus s for 
n = 1, 2, and 3. The dashed and dotted curves are the gas 
exponent for collapse with and without collisionless matter, 
respectively. The curves are computed according to tables 1 
and 2 with 7 = 5/3. For n = 2 the two curves overlap. The 
solid curve is the asymptotic exponent. <5, of the collisionless 
matter given by equations ([40Wl2) (Fillmore & Goldreich 
1984). 
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Figure 6. The shock location A s as a function of s. Figure 7. The trajectory and fluid variables of a particle 

as a function of time in spherical collapse with 7 = 5/3 for 
various s. The distance, r, and density, p, have been scaled 
by their respective values, r' ta , and, p' ta , at the turnaround 
time tt a , while the pressure, p, by its value, p s , immediately 
after the particle crossed the shock. 
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